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Abstract — Multiple-Input  Multiple-Output  (MIMO)  commu¬ 
nications  have  shown  great  promise  in  providing  high  spectral 
efficiency  for  wireless  ad  hoc  networks.  In  this  paper  we  study 
the  problem  of  joint  routing,  scheduling,  power  control  and 
bit  rate  selection  in  MIMO-based  ad  hoc  network  with  the 
goal  of  maximizing  the  system  throughput  that  satisfies  a 
given  end-to-end  traffic  demand.  We  formulate  this  cross-layer 
optimization  problem  as  a  linear  programming  problem.  It  is 
well  known  that  the  scheduling  problem  arising  in  wireless 
network  is  NP-hard  due  to  the  exponential  optimization  space. 
Accommodating  MIMO  links  in  this  work  greatly  increases  the 
computation  complexity.  The  cross-layer  optimization  problem 
can  be  solved  by  applying  an  iterative  approach  proposed 
in  [10],  and  the  original  problem  is  decomposed  into  sub¬ 
problems  and  solve  them  iteratively.  Numerical  results  are  given 
to  validate  the  solution. 

I.  Introduction 

Multiple-Input  Multiple-Output  (MIMO)  systems  have 
shown  great  promise  in  providing  high  spectral  efficiency 
for  single  user  wireless  link  without  interference  [1],  [2]. 
There  has  also  much  work  on  the  MIMO-based  cellular 
networks,  which  include  MIMO  multiple  access  (MIMO- 
MAC)  [3],  [4]  and  MIMO  broadcast  systems  [5],  [6].  Both 
systems  have  one  common  end  of  the  communication  link  - 
either  the  transmitter  of  MIMO-BS  or  the  receiver  of  MIMO- 
MAC.  There  has  been  great  interest  in  extending  the  MIMO 
communication  to  the  multi-user  systems  with  interference. 
The  transmission  scheme  of  each  user  depends  on  that  of 
other  users  since  the  interferences  at  each  user  depend  on 
all  the  transmit  covariance  matrices  in  the  network. 

A  Gradient  Projection  (GP)  algorithm  [7]  and  a  Quasi- 
Newton  (QN)  method  [8]  were  presented  to  provide  subopti- 
mal  solutions  subject  to  the  constant  power  constraint.  In  [9], 
we  propose  an  Integer  Programming  (IP)  method  to  study 
the  total  system  capacity  problem  for  the  MIMO  system  with 
interference.  Both  transmitters  and  receivers  are  assumed  to 
have  perfect  knowledge  of  channel  state  information  (CSI). 
Each  MIMO  link  is  decomposed  into  parallel  independent 
eigenchannel.  Then,  each  eigenchannel  is  represented  by  a 
set  of  logical  links,  and  each  logical  link  is  determined  by 
different  discrete  powers  and  discrete  data  rates.  Also,  the 
transmitters  of  each  MIMO  link  are  subject  to  the  total  power 
constraint.  The  solution  specifies  the  set  of  logical  links  that 
can  transmit  simultaneously. 

The  broadcast  nature  of  wireless  radio  causes  interference 
among  the  nodes  sharing  the  common  communication  me¬ 
dia.  It  is  well  known  that  the  optimal  scheduling  problem  in 
wireless  network  is  NP-hard.  The  optimization  space  grows 


exponentially  with  the  number  of  links.  For  example,  if 
single  channel  is  used,  then  one  must  maximize  over  a 
polytope  with  2L  extreme  points,  where  L  is  the  number 
of  links.  More  specifically,  we  define  an  assignment  to  be  a 
specification  of  which  links  are  transmitting  and  which  links 
are  not  transmitting.  A  schedule  is  the  convex  sum  of  assign¬ 
ments.  The  main  challenge  to  the  throughput  maximization 
problem  has  been  that  the  space  of  all  assignments  is  too 
large.  In  order  to  solve  the  scheduling  problem,  the  space 
of  assignments  must  be  reduced.  In  [10],  the  authors  have 
developed  the  iterative  techniques  to  compute  the  optimal 
scheduling  for  the  wireless  mesh  network  even  when  co¬ 
channel  interference  is  present. 

Allowing  the  schedule  optimization  to  optimize  over 
MIMO  links  results  in  a  higher  dimension  problem.  How¬ 
ever,  the  main  results  and  the  algorithm  for  optimal  schedul¬ 
ing  from  [10]  are  still  available.  In  this  paper,  we  study 
the  problem  of  resource  allocation  in  MIMO-based  Ad  Hoc 
networks  by  jointly  determining  routing,  scheduling,  power 
control  and  bit  rate  selection  with  the  goal  of  maximizing 
the  system  throughput  that  satisfies  an  end-to-end  traffic 
demand.  After  solving  the  optimization  problem,  the  assign¬ 
ment  specifies  which  links/eigenchannels  are  transmitting, 
their  transmit  powers  and  bit  rates.  Routing  is  also  solved 
as  a  by-product.  During  searching  for  a  new  assignment,  the 
dimension  of  the  integer  programming  problem  significantly 
increases  by  considering  MIMO  links. 

The  remainder  of  this  work  is  organized  as  follows. 
Section  II  briefly  introduces  some  related  work.  In  Section 
III,  an  integer  programming  based  algorithm  [9]  is  presented 
to  maximize  the  total  system  capacity  for  MIMO  links  with 
interference.  In  Section  IV,  problem  definition  are  given, 
and  an  iterative  approach  to  solve  the  resource  allocation 
problem  is  presented.  The  approach  focuses  on  reducing 
the  size  of  the  space  of  assignments  and  finding  a  new 
assignment.  Then,  Section  V  shows  the  results  from  numer¬ 
ical  experiments.  Finally,  concluding  remarks  are  given  in 
Section  VI. 

The  notations  in  this  paper  are  as  follows.  The  boldface 
denotes  matrices  and  vectors.  For  a  matrix  A,  A  denotes 
the  conjugate-transpose.  | Aj  denotes  the  determinant,  and 
Tr( A)  denotes  the  trace  of  a  square  matrix.  AAO  means 
that  A  is  a  positive  semi-definite  matrix.  A  is  a  Hermitian 
matrix  if  A  =  At.  I  denotes  the  identity  matrix  with 
the  appropriate  dimension  from  the  context.  E[-]  denotes 
the  statistical  expectation.  C(x  x  y)  denotes  the  complex 
space  of  x  x  y  matrix.  min(x,y)  and  max(x,y)  denotes 
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the  minimum  and  the  maximum  of  two  real  numbers  x  and 
y,  respectively. 


II.  Related  Work 

There  is  increasing  interest  in  the  MIMO-based  ad-hoc 
network.  In  [1 1],  Blum  investigates  the  MIMO  capacity  with 
interference  where  single-user  detection  is  assumed  at  the 
receiver.  Without  CSI  at  the  transmitter,  Blum  shows  that  the 
optimum  is  either  the  optimum  interference-free  approach, 
which  puts  equal  power  into  each  antenna,  or  a  singular 
mode,  which  puts  all  powers  into  a  single  antenna,  if  the 
interference  is  either  sufficiently  weak  or  sufficiently  strong. 

In  [7],  Ye  and  Blum  study  the  asymptotic  behavior  of 
the  system  mutual  information  and  propose  a  Gradient 
Projection  (GP)  algorithm  to  maximize  the  total  system 
capacity  subject  to  the  constant  power  constraint.  In  [8],  a 
Quasi-Newton  (QN)  method  is  proposed  to  solve  the  ca¬ 
pacity  problem  by  approximating  the  inverse  of  the  Hessian 
matrix  instead  of  computing  the  real  one.  However,  only  the 
suboptimal  solutions  can  be  found  due  to  the  non-convex 
nature  of  the  optimization  problem. 

Utility  optimization  of  wired  networks  was  pioneered 
by  the  seminal  works  of  Kelly  [12],  Low  [13].  The  same 
framework  has  been  extended  to  ad  hoc  networks  in  [14], 
[15].  The  main  obstacle  facing  throughput  maximization  is 
that  the  size  over  which  the  optimization  is  performed  is 
exponential  with  the  number  of  links  [16] — [  1 8].  To  reduce 
the  optimization  space,  one  approach  is  to  arbitrarily  restrict 
the  focus  to  a  subset  of  all  assignments,  which  is  fol¬ 
lowed  in  [18].  However,  this  approach  reduces  the  resulting 
throughput.  Another  approach  is  to  develop  a  heuristic  to 
determine  the  subset  of  considered  assignments  [19],  [20]. 
The  performance  of  these  methods  typically  decreases  as 
the  number  of  links  increases  [19],  [20].  Finally,  one  can 
take  a  brute  force  approach,  and  consider  the  entire  set  of 
assignments.  This  approach  is  taken  in  [16],  where,  due  to 
computational  difficulties,  the  largest  network  that  could  be 
solved  has  only  15  links.  However,  the  authors  in  [10]  has 
developed  the  iterative  techniques  to  compute  the  optimal 
scheduling  for  the  wireless  mesh  network  with  more  than 
2,000  links. 

III.  Capacity  of  MIMO  Links  with  Interference 

Without  considering  routing  and  scheduling  for  end-to- 
end  flows,  we  first  consider  the  capacity  of  a  MIMO  inter¬ 
ference  system  with  L  links  where  each  node  is  equipped 
with  rit  transmit  antennas  and  nr  receive  antennas.  Suppose 
all  MIMO  nodes  communicate  in  the  same  channel,  and 
cochannel  interference  is  presented  when  more  than  one  link 
transmit  simultaneously.  The  details  of  this  section  can  be 
found  in  [9].  Note  that  we  discuss  random  channel  matrix 
in  this  section,  but  we  only  consider  deterministic  channel 
matrix  when  we  include  routing  and  scheduling. 

A.  MIMO  Links  with  Interferences 

We  use  Hj j  €  C(nr  x  nt)  to  represent  the  channel  gain 
matrix  from  the  transmit  antennas  of  link  j  to  the  receive 


antennas  of  link  i.  The  complex  base-band  signal  vector 
received  by  the  receiver  node  of  link  i  is  given  by 

y;  =  \  /'.II  .-x-  +  X  V TT',/X/  +  n- 

j=t 


where  xt  represents  the  normalized  transmitted  signal  of  link 
i.  Pi  denotes  the  signal-to-noise  ratio  (SNR)  of  link  i,  pi:) 
denotes  the  interference-to-noise  ratio  (INR)  of  link  i  due 
to  the  interference  from  link  j,  and  n,  is  the  noise  vector. 
The  entries  of  H, -j  and  n,  are  independent  and  identically 
distributed  complex  Gaussian  random  variables  with  zero 
mean  and  unit  variance. 

Let  Q,  =  E[yLix\ ]  be  the  covariance  matrix  for  the 
transmit  signal  vector  x,.  Tr(Qj)  specifies  the  transmit 
powers  of  transmit  antennas  of  link  i,  and  we  have  Q,M) 
and  Tr(Qj)  <  1  where  1  represents  the  normalized 
maximum  power.  The  interference-plus-noise  of  link  i  is 
jjti  v  77-7!  l..,x,  +  n,.  which  is  Gaussian  distributed 
with  covariance  matrix  R,  =  I  +  J2j= 1  j^i  Pi  ,  n-./Qii;  y 
where  I  =  E(n;nj).  The  ergodic  capacity  [2]  of  link  i  can 
be  written  as 


a 


=  E  (tog2 


■  jRj 


where  the  expectation  is  taken  over  the  distribution  of  H. 

The  optimization  problem  to  maximize  the  system  capac¬ 
ity  can  be  defined  as 

max  ’X  WjCj  (la) 

i= 1 

subject  to  TY(Qj)  <  1 
QiYO. 


where  Wi  are  pre-assigned  weights.  Due  to  the  nonconvexity 
of  this  problem,  only  a  suboptimal  solution  can  be  found. 


B.  Capacity  of  MIMO  Links  via  SVD 

Assume  the  MIMO  interference  system  with  a  rich  scatter¬ 
ing  environment.  Singular  Value  Decomposition  (SVD)  can 
be  used  to  obtain  the  independent  channels  for  a  MIMO  link. 
Consider  a  MIMO  link  with  nr  x  tit  channel  gain  matrix  H 
that  is  known  to  both  the  transmitter  and  the  receiver.  The 
channel  gain  matrix  H  is  decomposed  into 

H  =  UAV'i’, 

where  the  nr  x  nr  matrix  U  and  the  rit  x  tit  matrix  V  are 
unitary  matrices,  and  the  nr  x  nt  matrix  A  is  a  diagonal 
matrix  of  singular  values  {cy}  of  H.  Due  to  considering  a 
rich  scattering  environment,  the  matrix  H  is  full  rank  and 
Rank  (H)  =  M  =  min(nr,nt). 

The  channel  gain  matrix  H,;  from  the  transmitters  of 
eigenchannels  of  link  j  to  the  receivers  of  eigenchannels 
of  link  i  (i  f  j)  shown  in  Fig  1  can  be  written  as 

"  r  i;:h,(v..  (2) 

where  U,  and  Vy  are  from  the  SVD  of  matrices  H(I 
and  Hjj.  Then,  Hfj  denotes  the  channel  gain  from  the 
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where  the  expectation  is  taken  over  the  distribution  of  H, 
crj,  is  the  singular  value  of  the  Ith  eigenchannel  of  link  i, 
and  Pi,  is  the  transmit  power  at  the  Ith  eigenchannel  of  link 
i.  The  details  can  be  found  in  [9].  From  the  proof,  we  know 
that  equation  (3)  is  very  accurate  when  L  and  M  are  large. 

In  summary,  the  channel  covariance  coefficient  DXmVn 
from  the  transmitter  of  the  nth  eigenchannel  of  link  y  to  the 
receiver  of  the  mth  eigenchannel  of  link  x  can  be  represented 
as 

t2 


D, 


0 

HI 


if x^v- 


if  x  =  y  and  m  =  n 
if  x  =  y  and  m^n 


Fig.  1 .  The  channel  gain  matrix  Hy  from  the  transmitters  of  eigenchannels 
of  link  j  to  the  receivers  of  eigenchannels  of  link  i 

transmitter  of  the  kth  eigenchannel  of  link  j  to  the  receiver 
of  the  Ith  eigenchannel  of  link  i. 

The  aggregate  interference-plus-noise  at  the  receiver  of 
the  Ith  eigenchannel  of  link  i  has  the  covariance  £■  written 


0= fc=i 

where  Pjk  is  the  transmit  power  at  the  kth  eigenchannel  of 
link  j.  Then,  the  approximately  ergodic  capacity  of  link  i 
can  be  written  as 

'  ''  “  ''iP  H 


C.  Integer  Programming  Based  Algorithm 

Suppose  that  there  are  K  modulation/coding  schemes  and 
S  transmit  powers  available  for  the  mth  eigenchannel  of 
link  x,  then  associated  with  each  MIMO  link  x  is  the  set 
of  logical  links  xm>fe)S  with  1  <  m  <  M,  1  <  k  <  K 
and  1  <  s  <  S.  We  define  an  assignment  v  as  a  set  of 
logical  links  that  can  transmit  simultaneously.  Specifically, 
an  assignment  v  £  {0,  iyLxMxKxS ,  where  vXmkB  =  1 
implies  that  the  MIMO  link  x  is  transmitting  at  the  mth 
eigenchannel,  the  kth  bit-rate  and  the  sth  normalized  power. 


We  define  7  as  the  buffer  size  to  reduce  the  sensitivity  to 
interference.  For  a  particular  px,  the  maximum  data  rate  of 
the  mth  eigenchannel  of  link  x  is 

C”ax  =  k>g2  + 

where  Pmax  =  1  is  the  normalized  maximum  power  and 
7  =  2 dB  is  used  in  this  work. 

For  simplicity,  we  select  I\  linearly  separated  data  rates 
denoted  by  [C™*/K  ...  C™*].  Then,  the  SINR  threshold 
to  achieve  the  kth  data  rate  over  the  mth  eigenchannel  of 
link  x  can  be  represented  as 

SINRf  =  2kxC™x/K  -  1. 

J'm,k 

Note  that  SINRtkm  depends  on  the  link  x.  Similarly,  the 
set  of  powers  can  be  selected  as  S  linearly  separated  transmit 
powers  denoted  by  [l/S1  ...  1], 

Let  RXm  k  s  be  the  data  rate  over  the  logical  link  xmtk,s- 
The  aggregate  interference-plus-noise  at  the  receiver  of  the 
mth  eigenchannel  of  link  x  has  the  covariance  (x  which 
can  be  written  as 

L  M  K  S 

Cxm  =  PX'V  Pyn,t,qVyn,t,q  +  ll 

y—l^y^Lx  n= 1  t=  1  q—  1 

where  PVn  t  is  the  transmit  power  of  logical  link  yn,t,q- 
The  logical  link  xm^k,s  to  be  active  must  satisfy 

"  />'  "  !  >  SINRlhm  k , 

Xm 

which  is  equivalent  to 

PxPxm,XmPxm.h,s  —  PINRxm  fcClEm  >  00  X  (vXm.,k,a  ~  l)  ■ 

(4) 

In  practice,  oo  is  replaced  with  some  large  number.  Define 
0{n )  and  I[n )  as  the  sets  of  links  that  are  outgoing  from 
and  incoming  to  node  n.  For  simplicity,  we  assume  that  only 
one  link  can  transmit  if  multiple  links  share  one  common 
node. 

Problem  (1)  can  be  approximated  by  solving  the  following 
integer  programming  problem. 

L  m  k  s 

max 

x=l  m= 1  k— 1  s=l 

subject  to 

PxDXm,XmPXm^  -  siNRihm  k(Xm  >  rXm  k  s  (yXmtht,  -  1) 

for  each  logical  link  xm^,s  (5a) 

m  k  s 

m=  1  k—1  s=l 

X]  Vx  -  1  for  each  node  n 

{xEO(n);x£l(n)} 

K  S 


xmk  a  VXm.  k  a  ^  1  f°r  eaCh  X  (5b) 


(5c) 


EE 

k= 1  s— 1 


<  vx  for  each  eigenchannel  of  link  x 


(5d) 


VXm,k,B  ^  {b)  1}  : 
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where 

r  Xm,k„=SINR?mkx 

L  M  K  S  \ 

Px,v  Dxm,yn  X/  X/  +  1  J  ' 

y—l^y^x  n=  1  t— 1  q— 1  / 

Constraint  (5b)  shows  that  the  total  power  at  the  transmitters 
of  link  x  must  be  no  more  than  the  normalized  power  1. 
Constraint  (5c)  enforces  that  only  one  link  can  be  active 
when  multiple  links  share  one  common  node  n.  Constraint 
(5d)  enforces  that  no  logical  link  of  each  eigenchannel 
of  link  x  can  be  active  if  link  x  itself  is  inactive.  One 
shortcoming  of  the  IP  algorithm  is  the  computation  difficulty, 
and  some  approximation  algorithms  based  on  Maximum 
Weighted  Independent  Set  (MWIS)  will  be  developed  in  the 
future  research. 

IV.  Resource  Allocation 

A.  Scheduling 

We  define  an  assignment  v  as  a  set  of  logical  links  that 
can  transmit  simultaneously.  In  other  words,  the  assignment 
specifies  which  links  are  transmitting,  which  eigenchannels 
are  active,  their  bit-rates,  and  their  transmit  powers.  Specifi¬ 
cally,  an  assignment  v  €  {0,  l}Lx Mx  K  xS,  where  vXm  k  s  = 
1  implies  that  the  MIMO  link  x  is  transmitting  at  the  mth 
eigenchannel,  the  kth  bit-rate  and  the  sth  normalized  power. 
The  set  of  considered  assignments  is  denoted  by  V,  while  the 
set  of  all  assignments  is  denoted  by  V.  The  size  of  V  is  the 
main  challenge  facing  optimal  scheduling.  Thus,  typically, 
we  only  consider  a  subset  of  all  assignments,  i.e.,  V  ^V. 
We  drop  the  terms  physical  and  logical  unless  it  is  unclear 
from  context. 

A  schedule  is  defined  as  a  convex  combination  of  assign¬ 
ments  and  can  be  represented  by  a  set  {av}  where  av  >  0 
and  Ev  av  <  1  where  av  is  the  fraction  of  time  that  the 
assignment  v  is  used.  Letting  R(v,xm)  be  the  data  rate 
over  the  mth  eigenchannel  of  link  x  during  assignment  v, 
the  average  data  rate  over  link  x  provided  by  schedule  { av } 
is  Ev“vEmfl(v,Ira). 

B.  Routing 

In  this  paper,  the  routing  of  MIMO-based  ad  hoc  networks 
is  studied  in  the  multicommodity  flow  model.  Let  $  be  a 
set  of  concurrent  sessions,  and  each  session  <f>  represents  a 
source-destination  pair  in  the  network.  A  traffic  demand  d0 
for  session  </>  has  to  be  transmitted  from  the  source  node  <f>s 
to  the  destination  node  <f>d  via  a  collection  of  paths.  That 
is,  the  flows  transmit  between  the  source-destination  pairs 
through  multi-path  routing.  All  links  are  directed;  bidirec¬ 
tional  communication  between  two  nodes  is  represented  by 
two  directed  links. 

We  denote  the  data  rate  associated  with  the  <jfh  flow  on 
link  x  by  r£,  and  define  O(n)  and  I(n )  as  the  sets  of  links 
that  are  outgoing  from  and  incoming  to  node  n.  At  each 
node  n,  which  is  the  source  of  flow  (j>,  the  following  flow 


Algorithm  1  Computing  an  Optimal  Solution 
1:  Select  an  initial  set  of  assignments  V(0),  set  i  =  0. 

2:  Solve  (6)  for  V  =  V(z)  and  compute  fi(i)  and  A (z), 
the  Lagrange  multipliers  associated  with  constraints  (6a) 
and  (6b),  respectively. 

3:  Search  for  an  assignment  v+  that  solves  (5) 

4:  if  the  assignment  v+  £  V(z)  or 

Ex=i  Px(i)  Em=i  R(.v+,xk)  <  A(z)  then 

5:  stop,  the  optimal  schedule  has  been  found. 

6:  else 

7:  set  V  (z  +  1)  =  V  (z)  U  v+,  set  i  =  i  +  1,  and  go  to 

Step  2. 

8:  end  if 


conservation  law  is  satisfied: 

5Z  rt  =  d<t> ■ 

xEO(n ) 

If  the  node  n  is  an  intermediate  relay  node  for  the  (f'h 
flow,  the  following  balance  equation  holds: 

E  4-  E  4  = 

x£0{n)  x£l(n) 

If  the  node  n  is  the  destination  for  the  (f  h  flow,  we  have 

-  rt  = 

x£l(n) 

C.  Resource  Allocation  Algorithm 

Suppose  the  MIMO-based  Ad  Hoc  network  with  L  links 
and  each  node  is  equipped  with  M  antennas.  With  this 
notation,  the  optimization  problem  of  resource  allocation  can 
be  described  as 

max  P 

subject  to: 

n :  source  of  <j) 

xEO(n) 

r%  -  ri  =  0’  n:  relay  of 

xEO(n)  x£l(n ) 

—  ^2  r£  =  —Pdtf,,  tv.  destination  of  </> 

M 

rt —  E  ax  ,  R{WiXm)  for  each  link  x  (6a) 

{<££<£}  vGV  m—  1 

5>v<l  (6b) 

vGV 

0  <  av  for  each  v  €  V.  (6c) 

(6)  is  a  linear  programming  problem,  which  is  extensively 
studied  in  [21]. 

In  theory,  (6)  is  solvable.  However,  there  is  a  significant 
computational  challenge  in  that  if  V  is  the  set  of  all  assign¬ 
ments,  then  the  vector  a.  has  2LMKS  elements.  Thus,  the 
size  of  the  space  over  which  the  optimization  is  performed 
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must  be  reduced.  This  idea  of  considering  a  reduced  space 
was  considered  in  [18]  and  [19],  however,  the  space  was 
constructed  arbitrarily.  In  this  work  the  space  is  constructed 
so  that  the  throughput  found  by  optimizing  over  the  reduced 
space  is  the  same  throughput  found  by  optimizing  over  the 
entire  space. 

The  idea  of  Algorithm  1  is  as  follows.  Given  an  initial 
set  of  randomly  selected  assignments  V,  solve  the  throughput 
optimization  problem  over  the  set  V  and  obtain  the  Lagrange 
multiplier  for  each  link.  Then,  find  a  better  assignment  v+  by 
solving  problem  (5).  If  v+  exists,  add  it  to  the  set  V  and  go 
back  to  solve  the  optimization  problem  again.  Otherwise,  the 
optimal  solution  is  obtained.  The  optimality  and  convergence 
property  are  proved  in  [10]. 

Because  problem  (6)  over  a  small  set  of  assignments  is 
a  linear  programming  problem,  the  computation  complexity 
of  Algorithm  1  is  determined  by  finding  a  new  assignment 
which  is  an  integer  programming  problem.  Specifically,  each 
MIMO  link  is  represented  by  a  set  of  logical  links.  The 
dimension  of  problem  (5)  significantly  increases  with  the 
number  of  MIMO  links.  Therefore,  Algorithm  1  is  fit  for 
moderate  size  of  MIMO-based  Ad  Hoc  networks. 

V.  Numerical  Experiments 

In  this  section,  we  present  some  numerical  results  to 
investigate  the  resource  allocation  in  MIMO-based  Ad  Hoc 
networks.  We  consider  a  randomly  generated  topology  with 
10  nodes  uniformly  distributed  in  a  square  region  of  1500m  x 
1500m  and  the  minimum  distance  between  any  two  nodes  is 
greater  than  300m.  Each  node  is  equipped  with  two  antennas 
and  the  maximum  transmit  power  for  each  node  is  Pmax=0.1 
watts  (20  dBm).  The  numbers  of  data  rates  and  powers  for 
each  eigenchannel  are  8  and  4,  respectively.  Path  loss  is 
represented  by  the  popular  two-ray  ground  reflection  model. 
The  total  bandwidth  in  the  network  is  20  MHz.  There  are 
two  sessions  in  the  network:  nodel  to  node  2,  and  node  9  to 
node  7.  For  simplicity,  the  traffic  demands  for  these  sessions 
are  [1  1]  Mbps.  The  computations  below  were  performed 
on  a  2.4MHz  AMD  FX-53  processor  with  8GB  RAM,  and 
CPLEX  vlO  is  used  to  solve  the  LP  and  ILP  problems. 

A.  Number  of  Iterations  until  Algorithm  1  Stops 

Figure  2  shows  how  the  throughput  increases  as  the  more 
assignments  are  added.  The  point  of  maximum  throughput 
occurs  when  the  stopping  condition  specified  in  Algorithm 
1  is  met.  Thus,  in  this  case,  Algorithm  1  stopped  after  43 
iterations. 

Note  that  only  one  assignment  is  added  at  each  iteration. 
Thus,  the  maximum  number  of  elements  in  V  is  the  number 
of  assignments  in  the  initial  set  of  assignments  plus  the 
number  of  iterations  required  by  Algorithm  1.  Hence,  we 
have  achieved  the  goal  of  determining  the  solution  to  (6)  for 
V  =  V  by  computing  the  solution  to  (6)  for  a  small  set  V. 
The  complexity  of  solving  linear  and  nonlinear  optimization 
problems  is  well  known,  and  is  not  investigated  here. 


Fig.  2.  Variation  in  the  computed  throughput  as  assignments  are  added. 


Flow  1 :  N1->N2 


Fig.  3.  Routes  and  flow  rates  for  Session  1. 


B.  Optimal  Routes  and  Flow  Rates 

Session  1  transmits  from  the  leftmost  node  to  the  right¬ 
most  node,  and  Session  2  transmits  from  the  topmost  node 
to  the  lowest  node.  The  routes  and  flow  rates  of  session  1 
and  2  are  shown  in  Fig.  3  and  Fig.  4,  respectively.  For  both 
sessions,  the  optimal  routes  are  multi-path  multi-hop,  and  the 
flow  rates  satisfy  the  flow  conservation  constraints.  Also,  it 
is  easy  to  see  that  bidirectional  communication  between  two 
nodes  is  represented  by  two  directed  links. 

VI.  Conclusion 

In  this  paper,  we  study  the  problem  of  resource  allocation 
in  MIMO-based  Ad  Hoc  networks  by  jointly  determining 
routing,  scheduling,  power  control  and  bit  rate  selection  with 
the  goal  of  maximizing  the  system  throughput  that  satis¬ 
fies  the  end-to-end  traffic  demand.  An  iterative  algorithm 
proposed  in  [10]  was  extended  to  solve  the  optimization 
problem  of  resource  allocation.  An  Integer  Programming 
based  algorithm  [9]  is  used  to  find  a  new  assignment  during 
each  iteration,  which  specifies  which  links/eigenchannels  are 
transmitting  and  their  transmit  powers  and  bit  rates.  After 
solving  the  optimization  problem,  joint  routing,  scheduling, 
power  control  and  bit  rates  in  MIMO-based  Ad  Hoc  net¬ 
works  are  founded  and  the  optimal  solution  is  obtained. 
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Flow  2:  N9->N7 


(m) 


Fig.  4.  Routes  and  flow  rates  for  Session  2. 

The  computation  complexity  of  the  Integer  Programming 
method  prevents  the  proposed  algorithm  working  for  large 
size  MIMO  networks.  In  the  future  work,  some  greedy 
approximation  methods  will  be  developped  to  trade  off  the 
achieved  performance  and  the  computation  complexity. 
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